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[i]  We  use  theory  and  global  ray  modeling  to  investigate  how  the  potential  of  gravity 
waves  to  transport  momentum  flux  globally  from  the  lower  atmosphere  into  the 
mesosphere  and  lower  thermosphere  (MLT)  varies  with  horizontal  wavelength  and 
ground-based  phase  speed.  Ray  modeling  is  performed  using  the  Gravity  Wave 
Regional  or  Global  Ray  Tracer  (GROGRAT)  interfaced  to  realistic  three-dimensional 
global  winds  and  temperatures  from  0  to  100  km  altitude,  specified  by  fusing  analysis 
fields  at  lower  altitudes  to  GCM  results  higher  up.  We  focus  on  gravity  waves  in  the  short 
10-  to  50-km  horizontal  wavelength  range  that  are  unresolved  by  global  models  and, 
according  to  theory,  can  transport  appreciable  momentum  flux  into  the  MLT.  Ray  results 
for  different  seasons  reproduce  some  of  the  limits  derived  from  simple  wave  theory:  that 
horizontal  wavelengths  shorter  than  10  1cm  tend  to  be  removed  by  vertical  reflection  or 
evanescence  at  the  source  and  slower  phase  speeds  are  more  prone  to  critical  level 
removal,  leading  to  a  preference  for  waves  with  longer  horizontal  wavelengths  and  faster 
ground-based  phase  speeds  to  reach  the  MLT.  These  findings  are  compared  to  the 
wavelength  scales  currently  resolved  by  satellite  limb  and  nadir  sounders,  highlighting 
wavelength  ranges  currently  measured  and  those  currently  unresolved.  A  road  map  is 
developed  for  how  current  and  future  satellite  measurements  can  be  combined  to  measure 
the  full  space-time  spectrum  of  gravity  waves  relevant  to  eddy  flux  deposition  and 
momentum  forcing  of  the  global  MLT.  In  particular,  recommendations  for  new  satellite 
measurement  strategies  that  fill  current  measurement  gaps  are  provided. 
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waves,  J.  Geophys.  Res.,  113,  D24104,  doi:10.1029/2007JD009682. 


1.  Introduction 

[2]  On  parameterizing  the  effects  of  breaking  gravity 
waves  (GWs)  in  a  global  numerical  model  for  the  first  time, 
Holton  [1982]  was  able  to  reproduce  the  wind  and  thermal 
structure  of  the  mesosphere  and  lower  thermosphere  (MLT). 
Since  then  it  has  become  established  that  GWs  are  the  main 
dynamical  driving  force  for  the  MLT  [McLandress,  1998; 
Fritts  and  Alexander,  2003],  However,  the  evidence  from 
modeling  sUtdies  is  indirect  and  based  on  matching  observed 
global  wind  fields  (e.g.,  zonal  mean  wind  climatologies) 
with  output  from  general  circulation  models  (GCMs). 

[3]  It  is  obviously  preferable  to  derive  wave-induced 
global  mean-flow  driving  directly  from  measurements  of 
actual  gravity  wave  motions  in  the  atmosphere.  However, 
measurement  of  gravity  waves  from  space  has  been  highly 
problematic  until  comparatively  recently  (see  Table  1  and 
Wu  et  al.  [2006]  for  an  overview).  Certain  satellite  remote¬ 
sensing  techniques  resolve  only  portions  of  the  full  space¬ 
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time  gravity  wave  spectrum,  and  it  remains  unclear  which 
wavelength  scales  from  which  measurements  are  most 
relevant  to  momentum  flux  deposition  in  the  MLT.  Thus, 
our  goal  in  this  paper  is  first  to  estimate,  as  best  we  can  from 
modeling,  which  gravity  wave  wavelengths  are  capable  of 
transporting  and  depositing  momentum  flux  in  the  MLT, 
and  then  to  assess  the  status  of  current  satellite  remote 
sensors  in  resolving  these  flux-carrying  portions  of  the 
gravity  wave  spectrum.  Our  purpose  is  to  provide  theoret¬ 
ical  guidance  for  the  design  of  next-generation  satellite 
remote-sensing  instruments  and  integrated  missions  capable 
of  resolving  the  relevant  gravity  wave  wavelength  ranges 
responsible  for  driving  the  MLT.  In  other  words,  is  it 
possible  with  state-of-the-art  technology  to  measure  all  the 
relevant  scales  of  GWs  from  space  and,  if  so,  which  remote¬ 
sensing  techniques  have  to  be  combined  to  achieve  this 
aim? 

[4]  The  relevant  scales  to  observe  are  determined  by  how 
GW  momentum  flux  varies  as  a  function  of  horizontal  and 
vertical  wavelength,  but  owing  to  the  lack  of  global  obser¬ 
vations  it  is  unclear  which  part  of  the  gravity  wave  spectrum 
conveys  which  fraction  of  the  total  vertical  flux  of  horizon¬ 
tal  momentum  into  the  MLT.  In  particular,  GWs  of  almost 
any  horizontal  wavelength  could  provide  sufficient  drag  to 
explain  the  temperature  and  wind  structure  of  the  MLT,  if 
they  occur  sufficiently  frequently  (see  below).  Therefore  we 
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Table  1.  Overview  of  Satellite  Instruments  Used  for  GW  Studies 


Technique 

Instrument 

Full  Name 

Selected  References 

Infrared  limb  sounding 

LIMS 

Limb  Infrared  Monitoring  of  the  Stratosphere 

Fetzer  and  Gille  [1994] 

CRISTA 

Cryogenic  Infrared  Spectrometers  and 

Telescopes  for  the  Atmosphere 

Eckermann  and  Preusse  [1999] 
Preusse  et  al.  [2002] 

Em  et  al.  [2004] 

CLAES 

Cryogenic  Limb  Array  Etalon  Spectrometer 

Preusse  and  Em  [2005] 

SABER 

Sounding  of  the  Atmosphere  using 

Broadband  Emission  Radiometiy 

Preusse  et  al.  [2006] 

HIRDLS 

High  Resolution  Dynamics  Limb  Sounder 

Alexander  et  al.  [2008] 

GPS  radio  occultation 

GPS-MET 

Tsuda  et  al.  [2000] 

CHAMP 

CHAllenging  Minisatellite  Payload 

de  la  Torre  et  al.  [2006] 

Microwave  saturated  limb  radiances 

UARS/MLS 

Upper  Atmosphere  Research 

Satellite/Microwave  Limb  Sounder 

Wu  and  Waters  [1996] 
McLandress  et  al.  [2000] 

Jiang  et  al.  [2004] 

AURA/MLS 

EOS-AURA/Microwave  Limb  Sounder 

Wu  and  Eckermann  [2008] 

Infrared  sublimb  view 

MSX/MWIR 

Midcourse  Space  Experiment/Mid- Wavelength 
Infrared  Radiometer 

Dewan  et  al.  [1998] 

Infrared  nadir  view 

AIRS 

Atmospheric  Infrared  Sounder 

Alexander  and  Barnet  [2007] 

Microwave  nadir  view 

AMSU 

Advanced  Microwave  Sounding  Unit 

Wu  and  Zhang  [2004] 
Eckermann  et  al.  [2006] 

should  ideally  design  our  instruments  to  measure  all  scales 
of  GWs  which  potentially  can  propagate  from  the  forcing 
altitude  into  the  MLT.  Fortunately  there  are  physical  limits 
on  the  range  of  wavelength  scales  for  freely  propagating 
GWs. 

[5]  The  intrinsic  frequency  Co  of  a  GW  is  limited  by  the 
relation 

N>\cb\>\f\,  (1) 

where  N  is  the  buoyancy  frequency  and  /  is  the  Coriolis 
parameter.  The  frequency  limit  (1)  implies  horizontal 
wavelength  limits.  The  Coriolis  parameter  inhibits  long 
horizontal  wavelength  GWs  at  high  and  mid  latitudes  and 
we  therefore  observe  an  increase  of  the  average  horizontal 
wavelength  toward  low  latitudes  [. Alexander  et  al.,  2002; 
Em  et  al.,  2004;  Wang  et  al.,  2005].  Experience  shows  that 
these  long  horizontal  wavelength  GWs  in  the  tropics  can  be 
observed  from  space  (e.g.,  Em  et  al.  [2004]  and  section  5). 
More  challenging  to  observe  are  waves  close  to  the  short 
wavelength  limit.  If  the  intrinsic  frequency  Co  approaches  the 
buoyancy  frequency  N,  the  wave  is  reflected  or  is 
evanescent  at  its  source  (the  latter  hereafter  referred  to 
simply  as  an  evanescent  GW).  There  is  a  minimum 
horizontal  wavelength  for  those  waves  which  can  propagate 
upward.  In  section  2  we  review  previous  work  [e.g., 
Schoeberl,  1985;  Iwasaki  et  al.,  1989;  Marks  and 
Eckermann,  1995;  Kim  et  al,  2003]  on  how  a  short 
horizontal  wavelength  boundary  can  be  estimated  from  the 
dispersion  relation.  This  wavelength  limit  depends  on  the 
intrinsic  frequency  of  the  GWs  and  on  the  buoyancy 
frequency  and  therefore  varies  with  altitude,  location  and 
season.  For  intrinsic  phase  speeds  slower  than  150  m  s-1 
the  horizontal  wavelength  limit  varies  between  10  and 
50  km  given  typical  stratospheric  conditions  [Alexander 
and  Dunkerton,  1999], 

[6]  The  primary  focus  of  this  paper  is  therefore  to 
investigate  the  short  horizontal  wavelength  boundary  due 
to  wave  reflection  and  wave  evanescence  and  its  variations 
in  a  realistic  three-dimensional  atmosphere  and  for  different 
seasons.  In  this  paper  we  restrict  attention  to  those  GWs  that 


can  freely  propagate  upward.  Of  course  it  would  be  more 
interesting  to  investigate  how  the  various  limits  affect  the 
GW-induced  background  wind  acceleration.  Flowever,  this 
requires  knowledge  of  the  source  distribution  and  on  global 
scale  the  shape  of  the  source  spectrum  and  the  intemiittency 
of  various  wave  sources  is  not  well  known. 

[7]  In  particular,  it  is  impossible  to  infer  from  the  avail¬ 
able  evidence  whether  very  short  horizontal  wavelength 
waves  (horizontal  wavelengths  \h  smaller  ~100  km)  or 
mesoscale  and  long  horizontal  wavelength  GWs  (horizontal 
wavelengths  A/,  longer  ~100  km)  dominate  the  driving  of 
the  MLT.  For  the  longer  horizontal  wavelengths  global 
information  has  been  provided  by  satellites  and  GW  resolv¬ 
ing  global  models.  However,  uncertainties  are  still  large.  For 
the  short  horizontal  scales  we  have  evidence  in  particular 
from  observational  case  studies  and  regional  modeling. 
Here  the  primary  source  of  uncertainty  is  to  extrapolate 
results  from  these  local  events  to  a  global-scale  momentum 
flux  budget.  In  the  following  paragraphs  we  give  a  brief 
overview  of  the  effects  of  GWs  in  global  models  and  of  the 
evidence  we  have  for  short  and  long  horizontal  wavelength 
GWs,  respectively. 

[8]  The  interaction  of  a  GW  with  the  background  flow  is 
governed  primarily  by  its  horizontal  phase  speed  and 
depends  only  weakly  on  it’s  horizontal  wavelength  [e.g., 
Hines,  1997;  Warner  and  McIntyre,  2001;  Fritts  and 
Alexander,  2003],  For  instance,  the  Hines  [1997]  parameter¬ 
ization  integrates  over  a  horizontal  wave  number  spectrum  to 
yield  a  single  “characteristic”  horizontal  wavelength  for  the 
GW  spectrum.  In  GCM  studies  typical  values  range  from 
~60  km  (shortest  values  used  by  Manzini  and  McFarlane 
[1998]  and  Akmaev  [2001]),  and  125  km  [Charron  et  al., 
2002],  to  several  hundred  kilometers  ( Manzini  et  al.  [1997] 
and  longest  values  used  by  Manzini  and  McFarlane  [1998] 
and  Akmaev  [2001]).  This  shows  that  longer  horizontal 
wavelength  GWs  could  conceivably  provide  sufficient  mo¬ 
mentum  to  drive  the  MLT.  However,  since  a  large  number  of 
partly  compensating  tuning  parameters  is  used  in  the  Hines 
scheme,  these  choices  do  not  provide  any  compelling  evi¬ 
dence  for  preferential  GW  horizontal  wavelengths  [Akmaev, 
2001;  Fritts  and  Alexander,  2003]. 
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[9]  The  Warner  and  McIntyre  GW  parameterization 
scheme  [Warner  and  McIntyre,  2001]  implicitly  integrates 
over  a  wide  range  of  intrinsic  wave  frequencies  Co,  thus 
including  both  short  and  long  horizontal  wavelength  GWs. 
Em  et  al.  [2005]  isolated  the  longer  wavelength  portion 
(A/,  >  100  km)  of  the  assumed  semiempirical  horizontal 
wavelength  spectrum,  which  reduced  the  momentum  flux  by 
about  a  factor  of  two  [see  Em  et  al,  2005,  Figures  2c  and 
2d],  indicating  that  the  two  horizontal  wavelength  ranges 
might  be  of  equal  importance  for  mean-flow  driving. 

[10]  Hamilton  et  al.  [1999]  enhanced  the  spatial  resolu¬ 
tion  of  a  GCM  to  a  degree  that  a  realistic  middle  atmosphere 
was  reproduced  without  employing  GW  parameterizations. 
The  highest  horizontal  grid  resolution  was  30  km.  With 
increasing  computer  power  more  models  will  be  run  at  these 
kinds  of  high  spatial  resolutions.  However,  spectral  analysis 
of  output  from  these  models  reveals  that  the  wave  momen¬ 
tum  flux  spectrum  is  still  not  fully  resolved  indicating  an 
ongoing  need  for  GW  drag  parameterization  [see  Kim  et  al, 
2003,  section  6c].  The  realism  of  the  gravity  waves  sponta¬ 
neously  generated  and  explicitly  resolved  in  high-resolution 
global  models  is  unclear  [e.g.,  Horinouchi  et  al.,  2003]. 
Thus,  before  high-resolution  model  results  can  be  interpreted 
in  order  to  understand  GWs  in  the  real  atmosphere,  they  need 
to  be  validated  by  global  measurements,  as  must  the  GW 
drag  parameterizations  used  in  low-resolution  models.  In  all 
cases,  satellite  measurements  of  GWs  are  absolutely  critical 
for  validation. 

[11]  Satellite  and  long-duration  balloon  estimates  of  GW 
momentum  flux  [Em  et  al.,  2004,  2006;  Vincent  et  al, 
2007;  Alexander  et  al.,  2008]  reveal  GWs  carrying  consid¬ 
erable  momentum  flux.  However,  GW  momentum  flux 
estimates  from  current  generation  satellites  still  suffer  from 
relatively  coarse  spatial  sampling  [Em  et  al.,  2004]  and 
accordingly  have  too  large  uncertainties  in  their  momentum 
flux  estimates  for  calculating  reliable  vertical  gradients. 
Therefore  mean  flow  acceleration  was  not  estimated  from 
these  data.  Nevertheless,  the  satellite  and  long-duration 
balloon  observations  clearly  indicate  that  longer  horizontal 
wavelength  GWs  are  important  for  the  total  GW  momentum 
flux. 

[12]  Short  horizontal  wavelength  GWs  (A/,  <100  km) 
have  fast  vertical  group  velocities  and  therefore  can  rapidly 
propagate  from  the  troposphere  into  the  MLT.  In  addition, 
the  vertical  flux  of  horizontal  pseudo  momentum  density 
of  a  GW  is  proportional  to  its  horizontal  wave  number 
(i.e.,  inversely  proportional  to  A/,).  Thus,  where  short 
horizontal  wavelength  GWs  are  present  in  the  atmosphere, 
they  can  theoretically  carry  significant  wave  momentum. 
Accordingly,  Alexander  and  Dunkerton  [1999]  found  that 
wave  reflection  can  have  a  large  influence  on  the  GW 
induced  acceleration  of  background  winds  in  the  MLT,  if 
a  spectrum  that  favors  short  horizontal  wavelengths  is 
assumed. 

[13]  A  prominent  source  for  short  horizontal  wavelength 
GWs  is  deep  convection.  High-resolution  simulations  of 
convection  and  squall  lines  show  that  isolated  deep  con¬ 
vective  systems  force  predominantly  short  horizontal  wave¬ 
lengths  [Fovell  et  al.,  1992;  Piani  et  al.,  2000;  Lane  et  al., 
2001],  Circular  wave  fronts  like  the  ones  predicted  by 
convective  models  have  been  observed  in  satellite  sublimb 
data  [Dewan  et  al.,  1998]  as  well  as  in  MLT  airglow 


imagery  [Taylor  et  al.,  1987;  Taylor  and  Hapgood,  1988; 
Sentman  et  al,  2003;  Suzuki  et  al.,  2007a],  However,  other 
models  and  mechanisms  of  GW  generation  by  deep  con¬ 
vection  predict  the  generation  of  longer  horizontal  wave¬ 
length  GWs  [Salby  and  Garcia,  1987;  Pfister  et  al.,  1993] 
and  such  waves  have  also  been  observed  in  situ  [Pfister  et 
al.,  1993]  and  from  satellites  [Preusse  et  al,  2001;  Jiang  et 
al.,  2004;  Preusse  and  Em,  2005], 

[14]  A  second  prominent  source  of  GWs  is  wind  flow 
over  orography  exciting  mountain  waves.  Mountain  waves 
have  been  extensively  studied  over  decades  (see,  e.g.,  the 
early  work  by  Queney  [1948]  and  reference  papers  by  Fritts 
and  Alexander  [2003]  and  Kim  et  al.  [2003]).  Depending  on 
details  of  the  topography  all  horizontal  scales  from  a  few 
kilometers  to  several  hundred  kilometers  can  be  forced  and 
have  been  observed  in  the  upper  troposphere  and  strato¬ 
sphere  [e.g.,  Nastrom  and  Fritts,  1992;  Dornbrack  et  al., 
2002;  Eckermann  et  al.,  2007;  Alexander  and  Barnet, 
2007], 

[15]  The  most  convincing  evidence  for  short  horizontal 
wavelength  is  provided  from  the  combination  of  air  glow 
images  with  radar  measurements  of  the  background  wind. 
These  data  have  been  evaluated  for  momentum  flux  in  the 
MLT.  Peak  flux  magnitudes  per  unit  mass  can  be  as  high  as 
400  to  900  m2  s-2  [Smith  et  al.,  2006;  Fritts  et  al.,  2002] 
and  accelerations  may  reach  up  to  80  m  s_1  in  1  hour. 
However,  nightly  or  monthly  averages  are  much  smaller. 
The  stations  measuring  momentum  flux  report  nightly  or 
monthly  averages  (peak  values)  of:  Starfire  Optical  Range 
in  New  Mexico,  30  m2  s-2  (10-100  m2  s-2)  [Swenson  et 
al.,  1999;  Tang  et  al.,  2002];  Shigaraki  in  Japan,  5  m2  s-2 
(20  m2  s-2)  [Suzuki  et  al.,  2007b];  Hawaii,  2  m2  s-2 
(nightly  averages  0-5  m2  s~2)  [Tang  et  al.,  2005];  Halley 
at  Antarctica,  6-10  m2  s-2  [Espy  et  al.,  2004];  and  Rothera 
at  Antarctica,  10-40  m2  s_  [Espy  et  al.,  2006].  The  local 
enhancement  above  the  Antarctic  Peninsula  (Rothera)  is 
also  found  at  lower  altitudes  for  longer  horizontal  wave¬ 
lengths  in  satellite  data  [e.g.,  Fetzer  and  Gille,  1994; 
Eckermann  and  Preusse,  1999;  Jiang  et  al.,  2002;  Em  et 
al.,  2004],  while  Hawaii  is  located  at  a  local  minimum  of 
subtropical  convectively  generated  GW  activity  [e.g.,  Jiang 
et  al.,  2004;  Ern  et  al.,  2004].  Assuming  these  short 
horizontal  wavelength  GWs  completely  deposit  their 
momentum  in  a  20-km-wide  vertical  altitude  region,  at 
Halley,  Rothera  and  New  Mexico  they  would  contribute  a 
significant  fraction  of  the  total  wave  drag.  At  some  places 
(e.g.,  Rothera,  New  Mexico)  the  observed  momentum  of  the 
short  wavelength  GWs  alone  might  suffice  to  drive  the  MLT, 
if  these  values  were  representative  of  the  zonal  means. 
However,  the  examples  also  show  that  global  observations 
are  required  for  the  interpretation  of  these  few,  isolated 
measurements. 

[16]  The  evidence  for  short  horizontal  wavelength  GWs 
stems  from  local  observations  and  modeling.  Estimating 
global  momentum  flux  from  these  isolated  sources  naturally 
implies  large  uncertainties.  Can  global  GW-resolving  mod¬ 
eling  with  high  resolution  GCMs  now  or  in  the  near  future 
close  this  gap?  Lane  and  Knievel  [2005]  show  that  the  GW 
momentum  flux  over  deep  convection  only  converges  at  a 
model  horizontal  grid  resolution  of  500  m.  In  addition,  the 
types  of  GWs  forced  by  convection  in  global  models 
depends  to  a  large  degree  on  the  ways  in  which  deep 
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convection  is  parameterized.  Different  parameterizations 
can  yield  enormous  differences  in  the  nature  of  the  forced 
waves,  suggesting  considerable  unreality  in  explicitly 
resolved  GWs  in  global  models  [e.g.,  Horinouchi  et  al., 
2003],  Thus  reliable  global  modeling  from  this  source  is 
at  present  infeasible. 

[17]  The  horizontal  wavelength  distribution  of  those  GWs 
which  drive  the  MLT,  if  known,  would  be  highly  relevant 
for  GW  drag  parameterizations.  Though  GWs  of  the  same 
phase  speed  and  amplitude  but  with  different  horizontal 
wavelengths  break  in  a  broadly  similar  fashion  within 
similar  background  wind  conditions  (the  basis  of  several 
GW  parameterizations,  see  above),  they  follow  different 
propagation  paths  from  low  to  high  altitudes  and  carry  and 
deposit  different  momentum  fluxes. 

[is]  Answering  this  question  about  the  most  relevant 
horizontal  GW  scales  requires  global  observations  of  the 
whole  spectrum  of  those  GWs  with  the  potential  to 
transport  wave  momentum  from  the  troposphere  or  tropo- 
pause  into  the  MLT.  True  global  coverage  is  only  provided 
by  satellite  instalments  and  an  increasing  number  of 
different  instruments  and  techniques  are  now  exploited 
for  GW  studies.  An  overview  of  GW  satellite  observations 
is  given  by  Wit  et  al.  [2006],  Table  1  summarizes  different 
techniques,  the  respective  instruments  and  selected  scien¬ 
tific  citations. 

[19]  Could  we  observe  all  the  upward  propagating  GWs 
in  the  middle  atmosphere  if  we  were  to  combine  observa¬ 
tions  by  different  remote-sensing  techniques?  In  order  to 
answer  this  question,  we  first  need  to  know  which  GWs  are 
able  to  propagate  upward  into  the  MLT.  Gravity  waves  can 
have  very  small  spatial  scales  and  some  waves  lie  outside 
the  visibility  limit  of  any  current  satellite  instrument. 
However,  GWs  with  very  small  scales  cannot  propagate 
far  vertically  in  the  atmosphere  and  therefore  have  very 
limited  potential  for  coupling  the  upper  troposphere  with  the 
mesosphere. 

[20]  In  this  paper  we  derive  theoretical  lower  limits  for 
the  horizontal  and  vertical  wavelengths  of  those  gravity 
waves  which  can  propagate  upward  into  the  MLT.  Previous 
studies  have  estimated  these  limits  using  highly  idealized 
wave  equations  and  climatological  background  atmospheres 
[e.g.,  Schoeberl,  1985;  Iwasaki  et  ah,  1989;  Kim  et  al., 
2003].  Here  we  perform  a  more  comprehensive  and  realistic 
modeling  study  which  traces  gravity  waves  globally  from 
the  troposphere  to  MLT  using  detailed  ray  tracing  and 
realistic,  three-dimensional  specifications  of  the  background 
global  atmosphere  from  0  to  100  km  altitude  for  4  days  in 
2003,  representative  of  spring,  summer,  fall  and  winter 
conditions.  We  investigate  latitudinal  and  seasonal  varia¬ 
tions  of  the  horizontal  and  vertical  wavelength  boundaries. 
We  use  the  limits  resulting  from  this  modeling  to  develop  a 
road  map  for  the  optimal  combination  of  instruments  on  a 
GW-dedicated  satellite  mission. 

[21]  Section  2  summarizes  basic  GW  equations  and  out¬ 
lines  the  physical  laws  limiting  GW  propagation.  Section  3 
describes  the  model  experiments  employing  the  Gravity 
wave  Regional  Or  Global  RAy  Tracer  (GROGRAT)  [Marks 
and  Eckermann,  1995;  Eckermann  and  Marks,  1997]  and 
section  4  infers  from  the  model  results  which  fraction  of 
waves  reaches  80  km  and  100  km  altitude,  respectively,  the 
reasons  why  rays  cease  propagating  upward,  and  how  this 


varies  with  latitude  and  season.  Section  5  discusses  the 
prospects  of  gaining  global  information  from  satellite  remote 
sensors  on  the  different  scales  of  GWs  which  potentially  can 
convey  momentum  from  the  troposphere  into  the  MLT,  and 
section  6  summarizes  the  results. 

2.  Basic  Theory 

[22]  Evanescence  and  reflection  of  GWs  are  governed 
by  the  GW  dispersion  relation,  which  is  also  the  basis  of 
the  GROGRAT  ray  tracer.  A  fairly  general  form  of  the 
dispersion  relation  for  GWs  takes  the  form  [Marks  and 
Eckermann,  1995] 


-.2  J^Af+ZV  +  o2)  ^ 

k2  +  m2+a2  ’  (2) 

where  k /,  =  ( k 2  + 12)112  is  total  horizontal  wave  number  ( k  and 
/  are  the  zonal  and  meridional  components,  respectively),  m 
is  vertical  wave  number,  N  is  background  buoyancy 
frequency,  /  is  inertial  frequency,  and  a  =  where  Hp  is 
density  scale  height.  The  GW  intrinsic  frequency  is 


Co  =  10  —  ku  —  Iv  =  k/,c,  (3) 

where  the  intrinsic  horizontal  phase  speed  is 

c  =  c  —  U  cos{q J).  (4) 

[23]  Here  U  =  {it  +  v2)1/2  is  the  total  horizontal  wind 
speed  (u  and  v  are  the  zonal  and  meridional  components, 
respectively),  to  is  ground-based  frequency,  c  =  tolkh  is  the 
ground-based  horizontal  phase  speed,  and  4>  is  the  angle 
between  the  horizontal  wind  direction  and  the  horizontal 
wave  propagation  direction. 

[24]  In  the  limit  of  short  horizontal  wavelengths  and 
higher  intrinsic  frequencies,  and  assuming  N2  C$>  f2,  then 
(2)  simplifies  to 


N2kj 

+  m2  +  a2 


4=>  m2 


N2 


-(.kl  +  a2). 


(5) 


[25]  An  upward  propagating  wave  experiences  varying 
winds  so  that  c  changes  according  to  (4),  as  well  as 
vertically  varying  buoyancy  frequency  N.  A  gravity  wave 
can  be  reflected  at  a  turning  level  where  m 2  — >  0,  which 
from  (5)  occurs  at  a  critical  horizontal  wavelength 


^ h.crit  —  2tT 


Vn2  —  a2c2 


(6) 


[26]  We  see  from  (6)  that  the  very  smallest  values  of  this 
critical  wavelength  occur  when  N2  W-  a2c2,  so  we  approx¬ 
imate  the  absolute  short  horizontal  wavelength  limit  as 


\crit  =  2^.  (7) 

[27]  The  expression  on  the  right  of  (7)  is  simply  the 
hydrostatic  expression  for  the  vertical  wavelength.  For  a 
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typical  tropospheric  value  of  N  =  0.01  s_1,  (7)  tells  us  that 
the  shortest  horizontal  wavelength  is  ^-0.628c  kilometers, 
for  c  values  in  meters  per  second. 

[28]  The  maximum  vertical  flux  of  horizontal  pseudo 
momentum  density  that  can  be  carried  by  a  GW  depends 
on  both  the  horizontal  and  vertical  wavelength.  Following 
Em  et  al.  [2004]  the  vertical  flux  of  horizontal  pseudo 
momentum  density  Fp  of  a  monochromatic  GW  is  propor¬ 
tional  to  its  vertical  and  inversely  proportional  to  its 
horizontal  wavelength  as  well  as  proportional  to  the  square 
of  the  temperature  amplitude  T, 

fp*y/-  (8) 

[29]  Since  the  amplitude  saturation  limit  Tmnx  due  to 
static  instability  is  also  proportional  to  the  vertical  wave¬ 
length  [Eckermann  and  Preusse,  1999] 


the  maximum  momentum  flux  that  can  be  carried  by  a  wave 
is  proportional  to  the  third  power  of  the  vertical  wavelength 
[Preusse  et  al.,  2006], 

Ep,max  «  =  “T  •  (10) 

A/,  m5 

[30]  Thus,  (10)  tells  us  that,  for  waves  of  a  given  peak 
temperature  amplitude,  the  largest  momentum  fluxes  are 
carried  by  those  waves  with  short  horizontal  and  long 
vertical  wavelengths.  Of  course  source  and  saturation  pro¬ 
cesses  also  affect  the  distribution  of  temperature  variance 
among  waves  of  different  horizontal  and  vertical  wave¬ 
lengths  [e.g.,  Tsuda  et  al.,  1991;  Bacmeister  et  al.,  1996], 

3.  Ray  Modeling 

[31]  A  full  description  of  the  GROGRAT  model  is  given 
by  Marks  and  Eckermann  [1995]  and  Eckermann  and 
Marks  [1997].  Briefly,  the  model  uses  the  nonhydrostatic, 
rotational  GW  dispersion  relation  of  the  form  in  equation  (2). 
The  resulting  ray-tracing  equations  take  into  account 
horizontal  as  well  as  vertical  refraction  of  the  wave  vector 
due  to  gradients  of  the  background  wind  and  temperature 
fields.  Temporal  variations  of  the  atmospheric  background 
can  also  be  explicitly  taken  into  account,  but  are  neglected 
here  since  a  static  background  atmosphere  is  assumed.  A 
wave  is  traced  upward  as  long  as  GW  physics  allows  the 
free  propagation  of  the  wave  and  as  long  as  the  ray 
remains  in  the  domain  where  a  background  atmosphere 
is  provided.  A  detailed  discussion  of  the  different  physical 
and  technical  conditions  which  cause  GROGRAT  to  ter¬ 
minate  a  ray  trace  are  discussed  below. 

[32]  Strictly  speaking  GW  ray  tracing  may  be  only 
applied  where  the  WKB  approximation  is  satisfied.  There¬ 
fore  GROGRAT  normally  uses  threshold  criteria  to  termi¬ 
nate  rays  when  the  WKB  assumption  is  violated.  This  does 
not  mean,  however,  that  in  the  real  world  the  wave  cannot 


further  propagate  and  many  ray  models  trace  waves  through 
such  limits  [e.g.,  Schoeberl,  1985],  The  WKB  approxima¬ 
tion  states  that  the  background  atmosphere  changes  slowly 
compared  to  the  wavelength  of  the  waves.  This  is  more 
likely  violated  for  fast  phase-speed,  long  vertical  wave¬ 
length  waves.  Therefore  the  standard  WKB  termination  of 
the  rays  poses  a  problem  in  our  study.  If  we  include  the 
WKB-terminated  waves  in  our  statistics,  we  bias  our  results, 
because  we  introduce  a  termination  criterion  which  would 
not  prevent  real  waves  from  further  propagating,  and  if  we 
exclude  these  waves  completely  from  the  statistics,  we  bias 
our  results,  because  the  WKB  criterion  will  more  often  be 
violated  for  reflected  waves.  Since  in  this  study  wavelength 
statistics  and  not  the  details  of  the  ray  path  is  our  primary 
focus,  we  decided  to  disregard  the  WKB  thresholds  and 
continue  the  wave  propagation  past  WKB  limits. 

[33]  Wave  amplitudes  along  the  rays  are  calculated  using 
conservation  of  vertical  wave  action  flux  densities  and  are 
dissipated  using  scale-dependent  radiative  damping,  back¬ 
ground  turbulent  diffusion,  and  one  of  several  available 
wave  saturation  schemes,  which  all  tend  to  have  largest 
effects  on  short  vertical  wavelength  gravity  waves  near 
critical  levels.  To  speed  up  the  calculations,  waves  slowly 
propagating  toward  their  critical  level  singularity  are  iden¬ 
tified  and  removed  using  a  minimum  wave  amplitude 
criterion  and  a  minimum  vertical  group  velocity  criterion 
(see  Appendix  A  for  the  sensitivity  of  the  results  to  the 
minimum  amplitude  criterion).  Waves  that  reflect  vertically 
are  identified  using  the  normalized  parameter 
which,  when  significantly  greater  than  unity,  identifies 
breakdown  of  the  WKB  approximation,  which  clearly 
occurs  as  m 2  — >  0.  See  Marks  and  Eckermann  [1995] 
for  further  details. 

[34]  The  aim  of  the  modeling  study  is  a  global  investi¬ 
gation  of  GW  propagation  with  a  particular  focus  on  short 
horizontal  wavelength  GWs.  The  model  setup  follows  that 
described  by  Preusse  et  al.  [2006],  but  uses  shorter  hori¬ 
zontal  wavelengths  at  the  source.  Rays  were  launched  either 
at  5  km  or  20  km  altitude,  in  eight  different  horizontal 
propagation  directions,  every  5°  latitude  and  every  20° 
longitude.  The  rays  are  launched  for  each  spectral  compo¬ 
nent  to  be  investigated.  One  spectral  component  experiment 
(SCE)  is  defined  by  the  ground-based  phase  speed  c  and  the 
horizontal  wavelength  \h.  Six  phase  speeds  of  c  =  3,  6,  10, 
20,  30  and  50  m  s~ 1  and  horizontal  wavelengths  of  A/,  =  10, 
25,  and  50  km  are  launched  at  both  source  altitudes.  In  order 
to  enhance  the  resolution  in  horizontal  wavelength  close  to 
A h.crih  additional  rays  were  launched  with  \h  =  40  km  at  the 
5  km  launch  level  and  with  A*  =  1 8  km  at  the  20  km  launch 
level  for  all  phase  speeds.  We  have  also  performed  some 
sensitivity  tests  on  the  choice  of  launched  phase  speeds.  The 
results  of  the  sensitivity  tests  can  be  found  in  Appendix  B. 

[35]  The  background  atmosphere  (wind,  temperature  and 
geopotential  height  fields)  is  specified  from  ECMWF  anal¬ 
ysis  data  for  altitudes  below  55  km  and  output  from  a 
thermosphere-ionosphere-mesosphere-electrodynamics  gen¬ 
eral  circulation  model  (TIME-GCM)  [Roble  and  Ridley, 
1994]  experiment  for  altitudes  40-100  km.  ECMWF  reanal¬ 
ysis  data  are  used  in  numerous  transport  studies  and  capture 
the  synoptic-scale  features  of  the  troposphere  and  strato¬ 
sphere  well  [Borsche  et  al.  ,2007 ;  Em  et  al.,  2007].  In  order  to 
reproduce  the  actual  atmospheric  state,  the  TIME-GCM  run 
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Figure  1.  Zonal  mean  zonal  winds  for  15  January,  15  April,  15  July,  and  15  October  2003.  The  data  are 
zonal  averages  of  the  wind  fields  used  for  the  ray-tracing  experiments  and  composed  from  ECMWF  data 
below  55  km  and  a  TIME-GCM  experiment  above  40  km  with  a  smooth  transition  at  the  overlapping 
altitudes. 


was  nudged  at  30  km  altitudes  to  NCEP  reanalyses  and 
radiation-forced  migrating  tidal  components  at  the  lower 
boundary  are  provided  from  the  GSWM  tidal  model  [ Hagan 
et  al,  1995].  The  GCM  was  run  continuously  from  January 


2002  to  December  2004  in  support  of  the  TIMED  mission 
[Oberheide  et  al,  2006]. 

[36]  At  the  overlapping  altitudes  of  40-55  km  a  smooth 
transition  from  ECMWF  to  TIME-GCM  fields  is  achieved 


Launched  at  20  km  Launched  at  5  km 


horizontal  wavel.  [km]  horizontal  wavel.  [km]  horizontal  wavel.  [km]  horizontal  wavel.  [km 


0  20  40  60  80  100 

percentage  of  launched  rays 

Figure  2.  Percentage  of  rays  which  reach  altitudes  of  80  km  or  higher  (left)  for  20  km  launch  altitude 
and  (right)  for  5  km  launch  altitude.  Values  for  (top)  15  January  and  15  April  2003  and  (bottom)  15  July 
and  15  October  2003. 


6  of  16 


D24104 


PREUSSE  ET  AL.:  WHICH  GRAVITY  WAVES  CAN  REACH  THE  MLT? 


D24104 


a)  Reason  for  ray-stop  July  Launched  @  20km 
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Figure  3.  Reason  why  rays  were  terminated  by  GRO- 
GRAT:  (a)  results  for  rays  launched  at  20  km  altitude  and 
(b)  results  for  rays  launched  at  5  km  altitude;  background 
atmosphere  is  for  15  July  2003.  For  details,  see  text. 


Reason  for  ray-stop  July  Launched  @  5km 
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by  calculating  a  weighted  mean  of  the  two  data  sets,  where 
the  height-varying  weight  of  the  ECMWF  data  is  1  at  40  km 
altitude  and  drops  linearly  to  zero  at  55  km  altitude  and  vice 
versa  for  the  TIME-GCM  data.  Figure  1  shows  the  resulting 
zonal  mean  zonal  winds  for  15  January,  15  April,  15  July 
and  15  October  2003.  The  wind  fields  reproduce  salient 
climatological  features  and  thus  demonstrate  that  this  fused 
data  set  realistically  represents  important  atmospheric  struc¬ 
tures  such  as  the  polar  jets,  the  wind  reversal  above  the 
summer  mesosphere  and  the  quasi  biennial  oscillation 
(QBO)  in  the  tropical  stratosphere. 

[37]  ECMWF  and  TIME-GCM  data  are  merged  after 
interpolation  onto  a  common  grid  of  2.5°  latitude,  3.75° 
longitude  and  41  pressure  levels  corresponding  to  an 
approximate  pressure  altitude  spacing  of  2.5  km.  These 


atmospheric  fields  are  ingested  into  GROGRAT  and  inter¬ 
nally  reinterpolated  onto  a  regular  geometric  height  grid 
every  2.5  km  from  0.5  km  to  100.5  km.  Latitude  coverage  is 
from  85°S  to  85°N. 

4.  Results 

4.1.  Spectral  Dependencies 

[38]  Figure  2  shows  globally  averaged  percentages  of  all 
the  rays  which  reach  at  least  80  km  altitude  and  thus  have 
potential  to  carry  GW  momentum  from  low  altitudes  into 
the  MLT.  The  x  axis  gives  the  horizontal  wavelength  at 
launch  altitude,  the  y  axis  gives  the  ground-based  horizontal 
phase  speed. 

[39]  From  criterion  (7)  for  N  ~  0.01  s-1  (0.02  s-1)  we 
would  expect  that,  for  a  10  km  horizontal  wavelength, 
GWs  with  intrinsic  phase  speeds  faster  than  ~15  m  s-1 
(30  m  s-1)  would  evanesce  [cf.  Marks  and  Eckermann, 
1995].  Thus  any  vertical  shear  in  background  winds  will 
tend  to  reflect  these  waves  more  easily.  Consistent  with 
these  inferences,  Figure  2  reveals  a  low  percentage  of  free- 
propagating  10  km  horizontal  wavelength  GWs  at  80  km 
altitude  in  all  seasons.  For  horizontal  wavelength  longer 
than  10  km,  slow  ground-based  phase  speeds  (<10  m  s~*) 
are  less  likely  to  reach  the  MLT  than  faster  phase  speeds. 
Figure  2  reveals  that,  in  all  seasons,  the  GWs  with  the 
highest  probability  of  reaching  the  MLT  are  those  with 
long  horizontal  wavelengths  and  fast  ground-based  phase 
speeds. 

[40]  Figure  3  plots  probability  distributions  of  the  dif¬ 
ferent  criteria  used  by  GROGRAT  to  terminate  ray  integra¬ 
tions  in  these  experiments.  These  criteria  are:  ( 1)  evanescence 
(m2  <  0  at  the  source);  (2)  vertical  reflection  (m2  <  0  away 
from  source);  (3)  reached  upper  boundary  (z  =  100  km); 
(4)  wave  amplitude  completely  dissipated;  (5)  wave  stalled 
near  critical  level  (vertical  group  velocity  criterion;  see 
section  3);  (6)  wave  propagated  outside  longitude-latitude 
range.  Only  options  3  and  6  are  purely  practical  rather 
than  physical  criteria.  (Option  4  is  a  proxy  for  critical-level 
removal.) 

[41]  In  Figure  3a,  only  a  few  rays  are  evanescent  at  the 
source  (except  for  \h  =  10  km,  c  =  50  m  s-1),  but  many  of 
the  short  wavelength  rays  soon  become  external  through 
refraction  and  are  reflected  vertically.  On  the  other  hand,  the 
slow  waves  (low  phase  speeds)  are  mostly  dissipated  or 
stalled  near  critical  levels  and  are  therefore  not  apt  to  carry 
significant  momentum  flux  to  high  altitudes  either.  For 
the  fast  waves  (c  =  30  m  s^1  and  c  =  50  m  s-1)  with 
A/,  =  50  km  the  number  of  waves  reaching  the  upper 
boundary  at  100  km  is  significantly  lower  than  those 
reaching  80  km  (see  Figure  2).  This  difference  highlights 
significant  deposition  of  GW  momentum  into  the  MLT 
between  80  and  100  km  altitude. 

[42]  Equation  (2)  predicts  that  the  horizontal  wavelength 
A ]J  crit  delimiting  the  evanescent  from  the  internal  GW 
wavelength  regime  is  inversely  proportional  to  N.  This 
yields  more  evanescent  rays  at  the  lower  launch  altitudes 
where  N  is  smaller  as  shown  in  Figure  3b.  Comparing 
Figures  3a  and  3b  we  find  a  general  increase  of  dissipated 
and  stalled  waves  for  slow  phase-speed  (c  <  10  m  s"1) 
GWs.  The  number  of  dissipated  waves  increases,  because 
the  waves  are  launched  below  a  frequently  present  wind 
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Figure  4.  Altitude-latitude  cross  sections  of  the  percen¬ 
tage  of  rays  in  an  altitude-latitude  bin.  Shown  are  values 
for  (left)  April  and  (right)  July,  (a,  b)  Winds  from  Figure  1, 
and  GROGRAT  results  for  phase  speeds  of  (c,  d)  3,  (e,  f)  10, 
and  (g,  h)  50ms“',  respectively.  Launch  level  is  20  km,  and 
values  are  integrated  over  all  four  horizontal  wavelengths. 


reversal  between  the  troposphere  and  stratosphere  that  also 
induces  geographical  variations  discussed  in  section  4.2. 
Similar  effects  are  seen  in  global  model  parameterizations 
[Siskind  et  al,  2003]. 

[43]  Having  identified  the  processes  largely  responsible 
for  filtering  different  parts  of  the  wave  spectrum  by  the 
background  atmosphere,  we  can  better  interpret  the  seasonal 
dependences  observed  in  Figure  2.  The  main  difference 
between  the  different  months  is  the  different  wind  velocity 
in  the  middle  atmosphere.  Figure  1  shows,  from  the  ground 
to  100  km,  the  zonal  mean  of  the  zonal  winds  which  we 
have  used  for  our  ray-tracing  experiment.  At  equinoxes  the 
wind  speeds  are  generally  lower  and  therefore  reflection  of 
high  phase  speed,  short  horizontal  wavelength  waves  at 
turning  levels  is  less  likely.  For  high  phase-speed  GWs  the 


atmosphere  is  therefore  more  transparent  at  equinoxes  than 
at  the  solstices. 

[44]  On  the  other  hand,  slow  phase-speed  waves  (c  < 
10  m  s')  can  be  refracted  to  longer  vertical  wavelengths 
and  higher  vertical  group  velocities  if  they  propagate 
against  the  strong  solstitial  midlatitude  stratospheric  jets 
[cf.  Alexander,  1998;  Preusse  et  al,  2006;  Eckermann  et 
al.,  2007].  Figure  2  shows  that  they  are  more  likely  to 
reach  80  km  altitude  at  solstices  than  at  equinoxes. 
However,  only  for  July,  20  km  launch  altitude  and 
horizontal  wavelengths  of  50  km  the  likelihood  is  close 
to  50%.  In  most  other  cases  the  likelihood  for  the  slow 
waves  to  reach  80  km  altitude  is  smaller  than  30%  and 
frequently  even  below  10%. 

[45]  Comparing  rays  launched  in  the  troposphere  (at  5  km 
altitude)  with  those  launched  above  the  tropopause  (at 
20  km),  the  probability  for  the  rays  launched  at  the  lower 
level  to  reach  80  km  altitude  is  nearly  a  factor  of  2  lower 
than  the  probability  for  the  rays  launched  at  20  km.  One 
important  controlling  factor,  as  discussed  earlier,  is  the 
lower  buoyancy  frequency  in  the  troposphere,  which  affects 
predominantly  the  short  horizontal  wavelength  waves  and 
therefore  shifts  the  boundary  between  external,  nonpropa¬ 
gating  waves  and  internal,  propagating  waves  toward  longer 
horizontal  wavelengths.  In  addition,  the  wind  reversal 
between  tropospheric  westerlies  and  stratospheric  easterlies 
in  summer  adds  strongly  to  the  critical-level  filtering  of 
slow  phase-speed  waves. 

4.2.  Latitude-Height  Distribution 

[46]  Different  phase  speeds  exhibit  strongly  differing 
transparency  distributions  with  height  and  latitude  in 
Figure  4.  Zonal  mean  results  for  three  different  phase 
speeds  integrated  over  all  horizontal  wavelengths  are 
shown  for  April  and  July.  For  low  phase  speeds,  high 
wind  speeds  are  favorable  for  GW  propagation  since 
waves  propagating  against  the  background  winds  are 
refracted  to  longer  vertical  wavelengths  and  therefore  have 
higher  probability  of  reaching  the  MLT.  Consequently,  a 
high  percentage  of  the  c  =  3  m  s~*  rays  can  reach  high 
altitudes  in  the  winter  polar  vortex  of  the  southern  hemi¬ 
sphere  (80°S-40°S,  July)  and  in  the  subtropical  jet 
(10°N-40°N,  July).  The  former  result  is  consistent  with 
observations  showing  strong  penetration  of  quasi-stationary 
mountain  waves  from  the  southern  Andes  into  the  middle 
atmosphere  during  July  [Jiang  et  al.,  2002],  for  example. 
Both  features  are  less  pronounced  in  April  (equinox).  In  the 
tropics  wind  filtering  due  to  the  QBO  can  be  discerned  in  the 
results.  This  is  particularly  pronounced  in  April.  The  first 
strong  wind  gradient  occurs  between  20  and  30  km  and 
removes  about  half  of  the  waves,  the  wind  reversal  between 
30  and  40  km  filters  nearly  all  the  remaining  waves.  This 
mechanism  is  well  known,  applying  in  a  similar  way  to 
Kelvin  waves  and  other  tropical  wave  modes.  Current 
estimates  [Dunkerton,  1997]  attribute  the  forcing  of  the 
QBO  at  about  equal  parts  to  large-scale  tropical  waves  and 
smaller-scale  GWs. 

[47]  For  high  phase  speeds  wave  dissipation  is  less 
important,  but  waves  can  be  either  evanescent  at  the  source 
or  reflected  vertically  at  a  turning  level.  For  both  months  the 
number  of  high  phase-speed  waves  (c  =  50  m  s-1)  strongly 
decreases  in  the  wind  gradients  below  the  stratospheric 
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Figure  5.  Histogram  of  the  globally  integrated  distribution 
of  the  vertical  wavelength  at  various  altitudes  for  rays 
launched  at  20  km  in  July.  Also  shown  are  percentiles  for 
70  km  altitude.  For  details,  see  text. 

wind  maxima,  where  these  waves  are  reflected,  and  the 
number  of  rays  penetrating  to  the  MLT  peaks  at  the  tropical 
latitudes  where  wind  speeds  are  low  in  comparison  to  the 
mid  and  high  latitudes. 

4.3.  Vertical  Wavelengths 

[48]  The  histograms  in  Figure  5  plot  probability  density 
functions  (PDFs)  of  ray  vertical  wavelengths  at  various 
altitudes  for  the  20  km  altitude  launch  spectrum.  The 
distribution  is  a  global  integral  for  July.  Only  waves  which 
reach  at  least  80  km  altitude  are  shown.  The  PDFs  are 
expressed  as  percentages,  such  that  the  integral  over  each 
histogram  is  always  100%.  Owing  to  refraction  by  the 
background  winds  the  histogram  peaks  shift  back  and  forth 
in  vertical  wavelength  with  altitude.  The  relative  importance 
of  long  vertical  wavelength  waves  is  largest  at  70  km 
altitude  near  where  the  solsticial  midlatitude  jets  peak  in 
both  hemispheres  in  Figure  4b.  At  45  km  altitude  there  is  a 
hemispheric  asymmetry  between  a  strong  winter  polar 
vortex  and,  at  this  altitude,  much  weaker  winds  in  the 
summer  hemisphere.  This  could  be  responsible  for  the 
double-peak  structure  of  the  45  km  histogram  in  Figure  5. 
For  the  70  km  histogram  we  have  also  highlighted  in 
Figure  5  the  5th,  30th,  50th,  70th  and  95th  percentiles;  that 
is,  the  sum  of  the  bins  left  of  the  5%  line  is  5%  of  the  total 
number  of  rays  at  this  altitude.  At  70  km  altitude  almost 
all  waves  have  vertical  wavelengths  between  4  km  (5th 
percentile)  and  48  km  (95th  percentile). 

[49]  On  the  basis  of  the  percentile  values,  Figure  6 
systematically  surveys  the  altitude  and  seasonal  dependence 
of  the  vertical  wavelength  distribution.  The  results  for 
20  km  launch  altitude  are  considered,  which  provide  a 
conservative  estimate  for  those  waves  coupling  the  tropo- 
pause  with  the  MLT.  To  better  distinguish  the  curves, 
5th,  50th  and  95th  percentiles  are  given  as  solid  lines, 
30  and  70th  percentiles  are  given  as  dashed  lines.  For 
most  curves  the  vertical  wavelength  steadily  increases 
with  altitude,  reaches  a  maximum  at  around  70  km, 
decreases  above  70  km  and  increases  again  above  90  km. 
Values  are  larger  at  the  solstices  than  at  the  equinoxes  owing 
to  stronger  solsticial  stratospheric  jets  that  refract  waves  to 
longer  vertical  wavelengths.  This  effect  would  be  stronger  if 
only  regions  of  strong  winds  (e.g.,  midlatitudes)  would  be 


considered.  The  seasonal  differences  are  most  noticeable  for 
the  large  percentile  values,  i.e.,  for  the  long  vertical  wave¬ 
lengths.  From  Figure  6  one  might  conclude  that  an  instru¬ 
ment  measuring  all  waves  between  3  km  and  35  km  vertical 
wavelength  captures  90%  of  the  waves  propagating  into  the 
MLT.  It  should  be  noted,  however,  that  there  is  a  certain 
arbitrariness  in  the  choice  of  wave  properties  at  the  source. 
While  low  vertical  wavelengths  are  removed  by  filtering 
processes  and  the  deduced  limit  is  therefore  well  represented 
by  the  chosen  launch  distribution,  the  inclusion  of  larger 
phase-speed  waves  would  shift  the  curves  toward  slightly 
longer  wavelengths. 

4.4.  Horizontal  Wavelength  Limits 

[50]  The  primary  aim  of  this  ray  modeling  is  to  quantify 
the  horizontal  wavelength  limits  imposed  on  short  horizon¬ 
tal  wavelength  GWs  by  evanescence  and  vertical  reflection 
and  to  compare  this  limits  to  the  sensitivity  ranges  of 
different  instrument  types.  From  Figure  2  we  have  seen 
that  waves  of  10  km  horizontal  wavelengths  very  rarely 
reach  the  MLT  and  that  fast  waves  of  50  km  horizontal 
wavelength  have  a  high  probability  of  reaching  80  km 
altitude.  Figure  2  shows  the  dependence  on  the  ground 
based  phase  speed.  However,  satellite  instruments  measure 
the  vertical  wavelength  of  the  waves  which  is  determined  by 
the  intrinsic  phase  speed.  In  the  following  paragraphs 
we  therefore  deduce  a  short  wavelength  boundary  that 
depends  on  the  vertical  wavelength  of  the  waves. 

[51]  Since  we  focus  on  the  limits  imposed  by  wave 
reflection,  we  have  chosen  to  launch  GWs  for  which  Xh  is 
of  similar  size  to  Xz.  For  larger  horizontal  wavelengths,  kh  = 
(k2  +  l2)1'2  becomes  much  smaller  than  \m\  and  can  be 
neglected  in  the  denominator  of  the  dispersion  relation 
(equation  (2)).  In  this  midfrequency  approximation  (i.e., 
f2  < C  w2  <C  N2),  the  vertical  wavelength  becomes  directly 
proportional  to  the  intrinsic  phase  speed  c  and  independent 
of  the  horizontal  wavelength.  For  a  50  km  horizontal 
wavelength  GWs  approximately  satisfy  the  midfrequency 
approximation  for  all  launched  phase  speeds  and  the  waves 


5.0  %,  50.0  %,  50.0  %,  70.0  %,  95.0  %  Percentiles 


verticol  wavelength  [km] 


Figure  6.  Altitude  profiles  of  the  5th,  30th,  50th,  70th,  and 
95th  percentiles  (see  Figure  5)  of  the  globally  integrated 
vertical  wavelength  distributions  for  the  four  different 
months. 
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Figure  7.  Histogram  of  the  globally  integrated  horizontal 
wavelength  distribution  (a)  at  70  km  and  (b)  for  altitudes 
between  25  and  50  km  and  vertical  wavelengths  between  12 
and  15  km  (blue  curves).  The  wavelengths  stay  close  to 
their  initial  values  and  the  percentiles  are  close  to  the 
wavelength  peaks.  Figure  7b  shows  results  of  the  rebinning 
to  only  four  bins,  the  fit  of  a  smooth  function  through  the 
initial  wavelengths  (red  crosses),  and  the  deduction  of  a 
border  wavelength  (vertical  dashed  line)  based  on  this  fit. 

with  A/,  =  50  km  are  therefore  representative  of  longer 
horizontal  wavelengths.  We  have  tested  the  dependence  on 
the  choice  of  horizontal  wavelengths  in  Appendix  B. 

[52]  Percentiles  are  a  simple  and  robust  method  for 
characterizing  the  onset,  center  and  width  of  smooth  dis¬ 
tributions.  Since  vertical  wavelengths  are  continuously 
refracted  by  vertically  varying  background  winds,  the 
resulting  histograms  (see  Figure  5)  are  quasi-continuous 
and  percentiles  can  be  applied.  The  method  does  not  work, 
however,  for  the  horizontal  wavelengths:  Horizontal  wave¬ 
lengths  are  refracted  by  horizontal  wind  shear.  However,  for 
the  short  horizontal  GWs  considered  in  this  study  the 
wavelength  at  higher  altitudes  remains  close  to  its  initial 
value.  This  is  shown  in  Figure  7a,  which  gives  the  distri¬ 
bution  of  horizontal  wavelengths  at  70  km  altitude  (blue 
curve).  The  percentile  values  stay  close  to  the  main  peaks, 
similar  to  initial  findings  at  65  km  using  the  COSPAR 
International  Reference  Atmosphere  (C1RA)  reported  by 
Eckermann  and  Marks  [1997].  The  noncontinuous  nature 
of  these  distributions  requires  a  different  analysis  method. 

[53]  The  method  we  use  to  infer  a  short  wavelength 
boundary  from  the  horizontal  wavelength  distributions 
consisting  of  several  sharp  peaks  is  demonstrated  in 
Figure  7b.  The  blue  curve  shows  a  horizontal  wavelength 
histogram  as  in  Figure  7a,  but  for  selected  vertical 


wavelengths  between  12  and  15  km  and  integrated  over 
altitudes  between  25  and  50  km.  The  histogram  shown  in 
blue  is  rebinned  to  four  coarse  wavelength  intervals 
around  the  initial  horizontal  wavelengths  and  the  maxi¬ 
mum  value  is  normalized  to  1.0  (black  rectangles).  Then 
the  function 


y  =  —  arctan(£x  +  rf)  +  -  (11) 

7 r  2 

is  fitted  to  the  data  using  the  free  parameters  (  and  r/.  The 
function  fits  a  smooth  transition  between  limiting  values  of 
0  and  1 .  The  dashed  vertical  line  shows  where  this  function 
(black  curve)  equals  0.4,  which  is  used  to  demark 
propagating  waves  from  evanescent  or  reflected  waves, 
hereafter  referred  to  as  the  boundary  wavelength  Xh  b.  (We 
have  also  fitted  a  parabolic  function  with  little  change  to  the 
resulting  \h  h  estimates). 

[54]  Figure  8  shows  how  Xh  b  varies  with  vertical  wave¬ 
length  in  the  25  to  50  km  altitude  interval.  Results  are 
shown  for  5  km  (dashed  dotted  lines)  and  20  km  (solid 
lines)  launch  altitude  for  the  four  months  (indicated  by 
color).  At  shorter  vertical  wavelengths  Xhb  is  also  short. 
This  agrees  with  equation  (2)  since  shorter  vertical  wave¬ 
lengths  correspond  to  lower  phase  speeds.  At  the  equinoxes 
the  atmosphere  allows  shorter  horizontal  wavelength  GWs 
to  propagate  upward  relative  to  the  solstices.  The  curves 
shift  to  longest  X/hb  values  in  July,  consistent  with  a  higher 
probability  of  wave  reflection  in  the  strong  solstitial  strato¬ 
spheric  wind  regimes.  It  should  be  noted  that  the  most 
prominent  MLT  structures  generally  explained  by  GW 
momentum  deposition,  such  as  the  cold  summer  mesopause 
and  the  zonal  wind  reversal  in  the  summer  MLT,  are  located 
above  these  strong  wind  regimes,  which  act  to  reflect  the 
fast  phase-speed  short  horizontal  wavelength  waves.  The 
largest  difference  is  between  a  launch  level  in  the  tropo¬ 
sphere  (5  km)  and  a  launch  level  above  the  tropopause 


Altitude  interval:  25-50  km 


horizontal  wavelength  [km] 


Figure  8.  Boundary  wavelength  Xh  b  between  waves 
vertically  propagating  from  the  launch  altitude  to  the  MLT 
and  evanescent  or  reflected  waves.  For  details,  see  text. 
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Figure  9.  Range  of  gravity  waves  propagating  from  low  to 
high  altitudes  delimited  by  the  red  and  purple  curves, 
compared  to  the  sensitivity  ranges  of  infrared  and  micro- 
wave  limb  sounding  and  nadir-viewing  satellite  instruments. 
The  pink  shaded  area  not  covered  by  any  of  these 
instrument  likely  can  be  partially  covered  by  a  future 
instrument  performing  sublimb  sounding  of  saturated 
radiance.  For  details,  see  text. 

(20  km),  which  is  also  consistent  with  our  earlier  dis¬ 
cussion  of  Figure  2. 


5.  Discussion 

[55]  Figure  9  compares  the  sensitivity  ranges  of  differ¬ 
ent  types  of  satellite  remote-sensing  instruments  to  the 
physical  limits  imposed  on  GW  horizontal  and  vertical 
wavelengths  as  determined  in  section  4.  From  Figure  8  we 
obtain  a  short  horizontal  wavelength  limit  \hb  of  ~20- 
25  km  for  vertical  wavelengths  greater  than  ~15  km  and 
\hb  ~  10-15  km  at  3  km  vertical  wavelength.  For  a 
short  vertical  wavelength  limit  we  use  the  5th  percentile 
from  Figure  6,  which  is  about  3  km.  Together  these 
values  form  the  red  line. 

[56]  The  intrinsic  frequency  of  GWs  is  limited  by  the 
Coriolis  parameter  such  that 


&  >f 


Cj 

1 


>  1. 


(12) 


[57]  This  does  not  yield  a  horizontal  wavelength  limit,  if 
the  ratio  \Cj/f\  can  approach  unity.  CRISTA  measurements 
indicate,  however,  that  the  ratio  \u>/f\  remains  larger  than 
1.4  in  general  [Preusse  et  al.,  2006],  at  least  for  vertical 
wavelengths  Az  >  5  km  which  CRISTA  can  resolve.  The 
purple  lines  in  Figure  9  show  the  corresponding  wave¬ 
length  limits  implied  by  a  limiting  value  of  \u>lf\  =  1.4  at 
25°  and  50°  latitude. 

[58]  The  red  and  purple  lines  in  Figure  9  therefore  delimit 
the  range  of  upward  propagating  GWs  which  have  the 
potential  to  transport  GW  momentum  from  the  troposphere 
or  tropopause  upward  into  the  MLT.  This  region  is  com¬ 
pared  to  the  sensitivity  ranges  of  the  three  remote-sensing 
techniques  most  frequently  used  to  observe  GWs  in  the 
middle  atmosphere. 


[59]  The  green  shading  in  Figure  9  gives  the  range  of 
wavelengths  resolved  by  limb  sounding  of  optically  thin 
emissions  (or  occultations)  based  on  the  analytic  radiative 
transfer  solution  derived  by  Preusse  et  al.  [2002], 


S(H,  k,m) 


1  dB  77 


4(72  +  a2 


(13) 


where  a  =  m/2REarth  =  27t/(2A :REartb),  ks  =  2tt/Xs  is  the  wave 
number  along  the  line  of  sight,  7  =  1/(2 HREartll),  REarti,  the 
Earth’s  radius,  H  the  pressure  scale  height,  and  is  the 
linear  expansion  in  temperature  T  of  the  blackbody  source 
function  B.  Temperature  retrievals  are  simulated  by  normal¬ 
izing  to  an  infinite  horizontal  wavelength  (S(IJ,  k,  m)/S(H, 
0,  m)).  This  analytic  solution  is  in  good  agreement  with 
explicit  end-to-end  numerical  simulations  for  the  CRISTA 
temperature  retrievals  [Preusse  et  al.,  2002]  as  well  as  with 
sensitivity  estimates  for  temperature  measurements  by  GPS 
radio  occultations  [Lange  and  Jacobi,  2003].  The  solid  green 
line  in  Figure  9  gives  the  horizontal  wavelength,  where 
the  sensitivity  for  retrieved  temperatures  (i.e.,  S(H,  k,  in)/ 
S(H,  0,  in))  drops  to  0.3.  The  sensitivity  is  estimated  on 
the  basis  of  the  assumption  that  the  GW  phase  fronts 
propagate  perpendicular  to  the  instrument’s  line  of  sight. 
This  is  the  most  unfavorable  geometry  for  resolving  the 
wave  and  so  it  seems  realistic  that  in  many  cases  waves 
with  half  of  that  horizontal  wavelength  could  be  observed 
[cf.  Preusse  et  al.,  2002].  This  is  indicated  in  Figure  9  by 
the  lighter-green  shading  between  the  solid  and  the  dashed 
green  line. 

[60]  The  maximum  achievable  vertical  resolution  for  limb 
sounders  depends  on  the  field  of  view,  the  vertical  sampling 
rate,  and  the  signal-to-noise  level  of  the  instrument.  CRISTA 
attained  about  3  km  vertical  resolution  in  its  final  temperature 
retrievals  [Riese  et  al,  1999;  Preusse  et  al.,  2002],  whereas 
SABER  and  FllRDLS  have  somewhat  better  vertical  reso¬ 
lutions.  Of  all  limb  sounding  instruments,  GPS  data  have  the 
best  vertical  resolution  of  ~1  km  in  standard  retrievals 
[Borsche  et  al. ,  2007]  and  even  significantly  higher  resolu¬ 
tion  achieved  in  specialized  retrievals  [Hocke  et  al.,  2003], 

[61]  A  second  class  of  instruments  traditionally  used  for 
GW  analyses  employ  microwave  limb  sounding  of  saturated 
radiance  [Wu  and  Waters,  1996].  We  here  use  the  2D 
sensitivity  function  for  these  observations  derived  by 
McLandress  et  al.  [2000],  but  for  half  the  width  of  the 
UARS-MLS  field  of  view,  and  show  the  range  where  this 
sensitivity  function  is  larger  than  0.3.  This  reduced  field  of 
view  is  reached  by  modern  microwave  sounders  such  as 
Aura  MLS  [Wu  and  Eckermann,  2008].  The  GW  wave¬ 
length  regime  covered  by  microwave  instruments  is  shown 
in  yellow  in  Figure  9.  Again  we  can  assume  that  at 
favorable  viewing  directions  relative  to  the  GW  propagation 
direction,  a  factor  2  shorter  horizontal  wavelengths  can  be 
observed  (light  yellow). 

[62]  Most  recently,  cross-track  scanning  nadir-viewing 
instruments  have  also  been  employed  for  GW  studies  as  a 
result  of  improvements  in  their  spatial  resolution  and 
radiometric  precision  [Eckermann  and  Wu,  2006].  Their 
horizontal  resolution  is  limited  by  the  size  of  a  measurement 
footprint,  whereas  their  vertical  resolution  is  limited  by  the 
width  of  the  weighting  functions  of  the  thermal  emissions 
they  measure  [Eckermann  and  Wu,  2006;  Alexander  and 
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Barnet ,  2007].  In  general,  these  instruments  are  sensitive 
only  to  waves  with  vertical  wavelengths  of  ~15  km  and 
longer.  Nadir  viewers  are  indicated  by  blue  shading  in 
Figure  9.  It  should  be  noted  that  for  all  these  techniques 
the  wavelength  boundaries  shift  toward  shorter  wavelengths 
when  observing  very  large  amplitude  events,  since  these 
remain  visible  even  if  the  instrument’s  sensitivity  is  well 
below  the  assumed  threshold  of  0.3. 

[63]  Figure  9  reveals  that  together  these  three  instrument 
types  cover  a  reasonably  broad  range  of  the  wavelengths 
typical  of  the  vertically  propagating  GWs.  However,  there 
remains  an  important  observational  gap  for  GWs  with 
horizontal  wavelengths  <60  km  and  vertical  wavelengths 
<15  km  indicated  by  pink  shading  in  Figure  9.  Numerical 
modeling  indicates  that  convective  forcing  via  the 
mechanical  oscillator  [Fovell  et  al.,  1992;  Piani  et  al., 
2000]  preferentially  generates  GWs  in  this  range.  This  gap 
can  likely  be  fdled  by  sublimb  remote  sensing.  For  instance, 
Dewan  et  al.  [1998]  show  two  cases  of  concentric  wave 
structures  with  horizontal  wavelengths  of  25  and  50  km  and 
(presumably)  vertical  wavelengths  of  16  and  18  km,  respec¬ 
tively,  emitted  from  mesoscale  storm  systems  as  observed  by 
the  sublimb  MSX  instrument.  The  sensitivity  of  sublimb 
sounders  can  be  estimated  by  adapting  the  approach  of 
McLandress  et  al.  [2000]  for  estimating  the  MLS  sensitivity 
but  applied  to  sublimb  angles  between  5°  and  30°  and  a 
200  m  width  of  the  field  of  view  (lcr  width  assuming  a 
Gaussian  shape).  The  results  show  that  the  complete 
range  of  horizontal  wavelengths  between  10  km  and  the 
sensitivity  limit  of  limb  sounding  can  be  resolved. 

[64]  A  multi-instrument  monitoring  strategy  from  space 
for  all  the  free-propagating  GW  wavelength  scales  in  the 
middle  atmosphere  therefore  seems  to  be  within  reach. 
Future  instruments  should  focus  on  improved  sampling 
density  across  as  well  as  along  the  orbital  track.  High- 
resolution  along-track  CRISTA  sampling  permitted  absolute 
values  of  GW  momentum  flux  to  be  inferred  from  temper¬ 
ature  fluctuations  [Em  et  al.,  2004,  2006]  and  HIRDLS  data 
will  provide  improved  estimates  due  to  this  instrument’s 
shorter  along-track  horizontal  sampling  distance  [ Alexander 
et  al,  2008],  Furthermore,  vector  estimates  of  GW  momen¬ 
tum  flux  (i.e.,  estimates  of  the  momentum  flux  value  as  well 
as  wave  propagation  direction)  would  be  possible  using 
infrared  limb  imaging  instruments  [Riese  et  al.,  2005].  In 
addition,  radiative  transfer  simulations  show  that  the 
amplitudes  of  atmospheric  GWs  are  underestimated,  if 
the  retrievals  cannot  take  into  account  the  horizontal  fine 
structure  of  the  GWs  [Preusse  et  al.,  2002].  In  current 
estimates  of  absolute  values  of  momentum  flux  [Em  et  al., 
2004,  2006;  Alexander  et  al.,  2008]  this  is  one  of  the  leading 
error  sources  and  uncertainties  can  be  as  large  as  a  factor  of  2 
[Em  et  al.,  2004].  Only  if  the  instrument  is  viewing  along 
track  and  samples  the  waves  adequately  can  this  underesti¬ 
mation  be  avoided  [cf.  Em  et  al.,  2004,  section  4.3], 

[65]  Having  simultaneous  information  on  the  vertical  as 
well  as  the  horizontal  wave  structure  permits  GW  charac¬ 
terization  that  then  allows  back-tracing  studies  to  identify 
lower  atmospheric  wave  sources  [Hertzog  et  al.,  2001; 
Wrasse  et  al,  2006].  For  those  waves  in  the  wavelength 
regime  where  nadir  and  limb  sounding  sensitivities  overlap, 
this  allows  for  the  combination  of  the  horizontal  wave  image 
gained  by  the  nadir  viewer  and  the  vertical  wavelength 


and  phase  information  obtained  from  the  limb  viewing 
instrument,  which  yields  enhanced  three-dimensional  GW 
characterization  and  instrument  cross  validation. 

6.  Summary 

[66]  In  this  paper  we  have  studied  the  transparency  of  the 
atmosphere  to  gravity  waves  of  short  horizontal  wave¬ 
lengths  using  theory  and  ray  modeling.  These  waves,  where 
present  in  the  atmosphere,  have  a  high  potential  for  trans¬ 
porting  and  depositing  eddy  momentum  into  the  meso¬ 
sphere/lower  thermosphere  (MLT)  region,  since  GW 
momentum  flux  is  inversely  proportional  to  the  horizontal 
wavelength.  The  very  shortest  horizontal  wavelength 
waves,  however,  are  evanescent  at  the  source  or  are  verti¬ 
cally  reflected  by  wind  shear.  Independent  of  the  horizontal 
wavelength,  higher  phase-speed  GWs  attain  longer  vertical 
wavelengths  and  can  potentially  transport  larger  momentum 
fluxes,  since  the  momentum  flux  of  saturated  waves  is 
proportional  to  the  vertical  wavelength  cubed. 

[67]  In  this  paper  we  have  performed  global  ray-tracing 
experiments  through  a  representative  global  atmosphere  at 
four  specific  times  (15  January,  15  April,  15  July  and 
15  October  2003)  for  a  source  of  6  different  ground- 
based  horizontal  phase  speeds  and  four  different  horizontal 
wavelengths  located  at  two  different  launch  altitudes  of 
5  km  and  20  km.  This  allows  us  to  investigate  which  GWs  on 
average  reach  the  MLT,  which  waves  do  not  and  the  physical 
reasons  why.  While  waves  with  high  phase  speeds  and 
longer  horizontal  wavelengths  (A/,  >  50  km)  nearly  always 
reach  the  MLT,  waves  with  very  short  horizontal  wave¬ 
lengths  are  often  reflected  or  are  evanescent  at  their  source 
and  waves  with  low  phase  speeds  are  frequently  dissipated 
near  critical  levels. 

[68]  Wave  reflection  occurs  in  regions  with  high  wind 
velocities,  in  particular  in  the  winter  polar  and  the  summer 
subtropical  stratospheric  jets  at  solstices.  This  is  particularly 
interesting  since  the  cold  summer  mesopause  and  the  wind 
reversal  in  the  MLT  are  located  above  such  jets  where  the 
short  horizontal  wavelength  waves  are  filtered.  Wave  dissi¬ 
pation,  by  contrast,  is  found  in  weak  wind  regimes  and  in 
particular  in  the  vicinity  of  wind  reversals. 

[69]  A  horizontal  wavelength  boundary  \h  b  of  ~2- 
30  km  demarks  those  waves  with  a  high  likelihood  to 
propagate  upward  (A/,  >  \hth)  and  those  likely  to  be  evanes¬ 
cent  at  the  source  or  vertically  reflected  ( \h  <  A Aj6).  \/,j,  varies 
with  season  and  tends  to  have  higher  values  at  solstices  and 
for  low-altitude  launch  levels. 

[70]  No  single  satellite  instrument  can  resolve  the  whole 
wavelength  range  of  all  the  propagating  GWs  in  the  middle 
atmosphere,  but  a  synergetic  use  of  different  kinds  of 
satellite  observing  techniques,  such  as  optically  thin  and 
saturated  limb  observations  as  well  as  saturated  sublimb  and 
cross-track  nadir  emissions  can  approach  the  goal  of  global 
monitoring  the  complete  3D  wavelength  range  of  all  those 
GWs  potentially  carrying  momentum  from  the  troposphere/ 
tropopause  region  into  the  MLT. 

[71]  Largest  improvements  in  satellite  GW  research  in  the 
future  are  to  be  expected  from  improved  analysis  techniques 
and  synergetic  use  of  multi-instrument  data  from  many 
satellites,  or  even  preferentially,  from  a  dedicated  GW 
mission  combining  several  such  instruments  on  the  same 
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Figure  Al.  Comparison  between  (left)  1%  relative  ampli¬ 
tude  threshold  used  in  standard  runs  and  (right)  1000  times 
reduced  threshold.  For  details,  see  text. 

platform  for  synchronized  common-volume  measurements. 
High  priority  in  instrument  development  should  be  given  to 
increased  spatial  sampling  both  along  and  across  track.  This 
can  be  gained,  for  instance,  by  proposed  new  infrared  limb 
imagers  providing  a  fully  resolved  three-dimensional  image 
of  the  wave  fronts  from  a  single  instrument. 

Appendix  A:  Dependence  on  the  Dissipation 
Threshold 

[72]  When  a  gravity  wave  approaches  a  critical  level, 
i.e.,  when  the  background  wind  speed  approaches  the 
ground  based  phase  speed,  then  the  intrinsic  frequency, 
the  vertical  group  velocity,  the  vertical  wavelength  and 


10  20  30  40  50  10  20  30  40  50 


the  saturation  amplitude  all  approach  zero,  too.  Since  the 
vertical  group  velocity  becomes  small,  a  ray  tracer 
approaches  the  critical  level  in  very  fine  steps  causing 
large  computational  costs.  However,  very  close  to  the 
critical  level  the  effects  of  such  a  wave  are  very  small 
because  of  the  small  saturation  amplitude.  Therefore  in 
the  standard  runs  of  this  paper  waves  are  terminated 
when  the  actual  amplitude  of  the  wave  falls  below  1% 
of  the  launch  amplitude.  In  this  appendix  we  show  that 
this  affects  waves  close  to  critical  levels  and  does  not 
remove  rays  with  small  amplitudes  in  the  lower  atmo¬ 
sphere  which  could  grow  to  large  amplitudes  in  the  MLT 
causing  large  accelerations  there. 

[73]  GROGRAT  uses  two  termination  criteria  which 
compete  below  a  critical  level:  waves  that  fall  below  a 
minimum  relative  amplitude  with  respect  to  the  launch 
amplitude  are  called  “dissipated,”  waves  with  group  veloc¬ 
ities  below  a  threshold  are  called  “stalled”  (see  sections  3 
and  4.1).  In  Figure  Al  we  compare  an  experiment  with  a 
minimum  relative  amplitude  threshold  of  1%  (Figure  Al, 
left)  to  an  experiment  with  a  1000  times  smaller  threshold 
for  July  and  20  km  launch  level  (Figure  Al,  right).  Figure  Ala 
shows  the  percentage  of  all  rays  reaching  at  least  80  km 
altitude  (see  Figure  2),  Figures  Alb  and  Ale  show  the 
percentage  of  stalled  and  dissipated  waves,  respectively  (see 
Figure  3),  and  Figure  Aid  shows  the  sum  of  stalled  and 
dissipated  rays.  There  are  some  minor  changes  (one  color 
step  or  less)  in  the  number  of  waves  reaching  80  km  altitude 
or  more.  However,  the  relative  importance  of  the  two 
termination  criteria  shifts  for  the  run  with  reduced  minimum 
amplitude  threshold.  In  particular  for  the  short  horizontal 
wavelength  waves  the  number  of  dissipated  waves  is 
significantly  reduced,  but  this  is  compensated  by  an  equally 
increased  number  of  stalled  waves  (changes  in  the  sum  of 
both  termination  criteria  are  again  very  minor).  The  exper¬ 
iment  therefore  confirms  that  the  minimum  amplitude 
criterion  indeed  removes  waves  below  critical  levels  and 
does  not  filter  other  waves  of  low  amplitudes,  the  latter 
potentially  having  large  impact  on  the  MLT. 

Appendix  B:  Dependence  on  the  Chosen  Launch 
Distribution 

[74]  In  this  study  we  are  primarily  interested  in  the  effects 
of  wave  reflection  and  wave  evanescence  on  the  ability  of 


horizontal  wavelength  [km] 


Figure  Bl.  As  in  Figure  7b  but  including  mesoscale  GWs 
of  100  km  and  200  km  horizontal  wavelength  at  launch 
level. 
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Altitude  interval:  25-50  km 


Figure  B2.  Comparison  of  the  boundary  wavelength  Xbb 
between  (thin  lines)  the  original  experiment  and  (bold  lines) 
the  experiment  extended  by  two  mesoscale  components  of 
100  km  and  200  km  horizontal  wavelength. 

waves  to  propagate  upward  from  sources  in  the  troposphere 
or  at  the  tropopause  toward  the  MLT.  We  therefore  chose 
wave  components  for  which  horizontal  and  vertical  wave¬ 
length  are  of  comparable  size.  However,  for  the  deduction 
of  the  boundary  wavelength  Xh  h  delimiting  the  propagating 
from  the  evanescent  or  reflected  waves  we  have  assumed 
that  at  50  km  horizontal  wavelength  the  waves  would 
satisfy  the  mid  frequency  approximation  and  that  these 
waves  are  representative  of  longer  horizontal  wavelength 
GWs.  In  the  strong  solsticial  winds  this  assumption  is  not 
always  valid.  We  therefore  test,  in  this  appendix,  whether 
the  choice  of  launched  waves  and,  in  particular,  the  omis¬ 
sion  of  true  midfrequency  components  with  mesoscale 
horizontal  wavelengths  affects  the  determined  boundary 
wavelength  Xh  b. 

[75]  In  section  3  we  describe  the  setup  of  the  ray-tracing 
experiments.  In  particular,  we  have  chosen  six  phase  speeds 
and  four  horizontal  wavelength  components.  We  now  extend 
these  experiments  by  two  further  horizontal  wavelengths  of 
Xb=  100  km  and  Xh  =  200  km.  Figure  B 1  is  for  the  same  data 
as  Figure  7b  supplemented  with  the  two  additional  horizon¬ 
tal  wavelengths  (see  section  4.4  for  a  detailed  description). 
Since  the  integral  of  the  histogram  (blue  curve)  is  normalized 
to  100%,  the  percentage  value  for  the  four  shortest  wave¬ 
lengths  is  now  smaller  than  in  Figure  7b.  The  maximum 
value  of  the  coarse  histogram  (shown  in  black)  is  normalized 
to  one  (y  axis  on  the  right)  and  therefore  Figure  7b  and 
Figure  B1  are  comparable.  The  arcus  tangens  curve  given  in 
grey  is  the  original  fit  from  only  four  horizontal  wavelength 
components  (as  in  Figure  7),  the  black  curve  shows  the  fit 
through  all  six  components.  It  can  be  seen  that  the  differ¬ 
ences  are  small. 

[76]  In  Figure  B2  we  compare  the  deduced  boundary 
wavelengths  X/hb.  The  thin  lines  reproduce  the  results  of  the 
fits  from  four  horizontal  wavelengths  shown  in  Figure  8. 
Overplotted  (bold  lines)  are  the  results  from  six  wavelength 
components.  For  most  vertical  wavelength  ranges  the  results 
are  almost  indiscernible,  in  particular  for  the  20  km  launch 
level.  For  the  5  km  launch  level  and  vertical  wavelengths 


around  15  km  the  horizontal  wavelengths  are  somewhat 
larger.  The  intraseasonal  variation  and  the  general  shape 
remains  unchanged.  However,  the  fits  for  six  horizontal 
wavelength  components  do  not  work  for  either  very  short  or 
very  long  vertical  wavelengths.  At  slow  phase  speeds  most 
waves  are  dissipated  or  stall  below  critical  levels.  Short 
horizontal  wavelengths  have  higher  group  velocities  and 
their  vertical  wavelengths  are  somewhat  longer  (see 
equation  (5))  and  therefore  the  midfrequency  waves  are 
even  more  susceptible  to  dissipation  than  the  short 
horizontal  wavelength  waves.  Therefore  the  number  of 
waves  decreases  at  longer  horizontal  wavelengths  and  the 
arcus  tangens  fit  fails.  Likewise,  only  very  strong  back¬ 
ground  winds  can  refract  the  two  longer  horizontal 
wavelength  components  to  very  long  vertical  wavelengths 
(larger  than  ~30  km)  whereas  shorter  horizontal  wave¬ 
length  waves  again  attain  longer  vertical  wavelengths 
according  to  (5).  Therefore  the  boundary  wavelength 
curves  determined  also  from  the  mesoscale  components 
cover  a  reduced  vertical  wavelength  range  but  otherwise 
confirm  the  results  we  have  obtained. 
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